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Abstract 

A  theoretical  analysis  of  the  scaling  law  and  the  value  of  the 
threshold  current  for  beam  plasma  discharge  (BPD)  is  presented,  based  on 
the  requirement  for  an  absolute  Instability  near  the  plasma  frequency. 

It  Is  shown  that  both  the  scaling  law  as  well  as  the  numerical  values  of 
Ic  are  consistent  with  the  experimental  data,  both  in  the  low  and  high 
pressure  regimes  for  weak  magnetic  field  experiments  (u>e  >  &e) .  The 
differences  in  scaling  between  the  two  regimes  are  attributed  to 
transition  from  Bohm  to  classical  diffusion.  It  is  found  that  the  value 
of  the  pressure  minimum  for  ignition  increases  linearly  with  the  ambient 
magnetic  field. 
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I.  Introduction 

Laboratory  studies  of  energetic  electron  beam  injection  experiments 
in  a  neutral  gas  filled  vacuum  chamber  carried  out  in  the  large  vacuum 
facility  at  the  Johnson  Space  Center  (JSC)  (Bernstein  et  al.  1978,  1979, 
1980)  have  provided  important  Informations  for  the  interpretation  of 
data  from  space  based  electron  beam  injection.  Perhaps  the  most 
important  aspect  was  the  determination  of  an  empirical  relationship  of 
the  form 


for  the  critical  current  Ic  required  for  beam  plasma  discharge  ignition 
and  the  values  of  beam  energy  (Eb),  ambient  magnetic  field  (B) ,  system 
length  L,  and  ambient  pressure  p.  The  value  of  \  •»  .5-1  and  the 
pressure  function  f(p)  was  a  function  with  a  minimum  at  pQ  «  20  uT  and 
varying  roughly  as  p**-’  to  p*^  above  and  below  this  pressure  (Kellog  et 
al.  1981).  The  minimum  value  of  Ic  *  10  mA  for  a  system 
with  L  «•  20  m,  -  1  keV  and  B  ■  1  G,  typical  of  the  experiment.  The 
operating  pressure  range  was  between  1-50  yT,  with  most  measurements  on 
the  1-20  uT  range.  Similar  results  have  been  reproduced  in  a  number  of 
other  experiments  (Konrad  et  al.  1982,  Lyakhov  et  al.  1982,  Bernstein  et 
al.  1983)  in  which  the  collislonal  ionization  by  the  beam  is  sufficient 
to  bring  the  plasma  density  of  the  system  to  the  point  that  <i>e  >  &e, 
where  <ue,  are  the  plasma  and  cyclotron  frequency.  The  observed  p,  B 
scaling  is  shown  graphically  in  Tig.  1. 

Triggering  of  BPD  has  been  long  associated  with  a  beam  plasma 
instability  between  the  electron  beam  and  the  beam  generated  plasma 
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(Kharchenko  et  al.  1962,  Getty  and  Smullln  1963,  Galeev  et  al,  1976, 
Llnson  and  Papadopoulos  1980,  Papadopoulos  1981,  Galeev  1983).  For 
o>e  >  and  finite  size  systems  Rowland  et  al.  (1981)  and  Papadopoulos 
(1982),  have  associated  the  triggering  of  the  BPD  with  the  threshold  for 
an  absolute  beam  plasma  instability  near  u)g.  As  explained  in 
Papadopoulos  (1981)  for  systems  such  as  the  JSC  tank  the  system  length 
(L  £  20  m)  is  not  long  enough  to  allow  convective  modes  to  grow  to 
sufficient  amplitude.  The  requirement  that  the  waves  grow  at 
frequencies  near  we,  is  connected  with  the  fast  rate  of  non-linear 
energy  transfer  of  the  beam  energy  to  ionizing  suprathermal  electron 
tails  (Papadopoulos  and  Coffey  1974,  Papadopoulos  1975,  Papadopoulos  and 
Rowland  1978,  Rowland  et  al.  1980,  Galeev  1983).  It  is  the  purpose  of 
this  paper  to  develop  a  model  for  BPD  ignition  and  the  expected  scaling 
laws  on  the  basis  of  the  criterion  for  an  absolute  instability  near 
u)e.  Notice  that  for  systems  with  u>e  <  0g  the  u)e  waves  are  in  the  lower 
hybrid  branch  (Manickam  et  al.  1975)  which  in  the  absence  of  internal 
wave  reflections  gives  always  convective  amplification.  Our  analysis 
therefore  applies  only  to  situations  where  u>e  > 

The  plan  of  the  paper  is  as  follows.  We  discuss  next  the  beam 
plasma  equilibrium  expected  on  the  basis  of  collislonal  ionization. 
Section  III  presents  the  instability  theory  for  the  configuration 
determined  in  section  II  and  derives  the  threshold  criteria.  Section  IV 
presents  a  comparison  of  the  model  to  the  BPD  ignition  values  determined 
in  the  JSC  experiment.  The  final  section  summarizes  the  findings  and 
discusses  their  applicability  to  other  situations.  The  numerical  values 
are  given  in  MKS  units  except  for  the  beam  energy  (keV),  pressure  (yT), 
magnetic  field  (G),  and  density  (cm“^). 
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II.  Pre-BPD  Penalty  Buildup 

Before  entering  the  instability  analysis  it  is  necessary  to 
establish  the  equilibrium  density  profiles  for  the  plasma  during  the 
collisional  ionization  stage.  We  assume  that  the  beam  density  profile 
is  given  by 


nb(r)  =  nbe 


2/  2 
-r  / a 


,-r2/a2 


eV.  wa 
D 


(i) 


where  1^  and  are  the  electron  beam  current  and  velocity  parallel  to 
the  magnetic  field  and  a  the  beam  radius  given  by 


V,  sin  0, 
d  a 
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(la) 


where  0^  is  the  equivalent  divergence  injection  angle  (Linson  and 
Papadopoulos  1980).  The  equation  for  the  ionization  at  midplane  (i.e. 
ignoring  the  z  dependence)  is 
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where  D  is  the  diffusion  coefficient,  L  the  system  length,  NQ  the 
ambient  neutral  density,  and  o  the  ionization  cross  section.  The  term 
^  n(r)  describes  the  axial  losses.  Notice  that  if  we  average  (2)  over 
the  volume  we  recover  the  zero  dimensional  description,  in  terms  of  the 
confinement  time  t  (Papadopoulos  1982),  i.e. 
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which  gives  the  steady  state  value  of  the  density  as 
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The  general  solution  of  eq.  (2)  In  terms  of  the  first  and  second  order 
Bessel  functions  IQ,  KQ  Is 


'  Vo  <s>  +  Vo  <s> 


(5) 


with 


h 
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(6) 


The  constants  Aj,  A2  to  be  found  from  the  boundary  conditions.  In  the 
thin  beam  limit  (a  «  b) ,  in  which  we  have  a  line  source  we  find 
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The  more  general  solution  gives,  for  regions  Inside  the  source 
(0  <  r  <  a) 
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while  outside  the  source  (a  <  r  <  b) 
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Guided  by  the  experimental  results  we  restrict  ourselves  here  to  the 

Vb  sin  9d  1/2  1/2 

thin  beam  limit  a  <  b,  which  Implies  that  - ^ -  <  - — -  .  In 

e 

this  case  the  input  paramets  to  the  instability  analysis  are  n^ ,  nQ,  a, 
and  b  given  by  eqs.  (1),  (6)  and  (7).  The  values  of  nb  and  nQ  In  the 
system  units  discussed  in  section  I  can  be  found  from  eqs.  (1)  and  (7) 
as 
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n 
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III .  Instability  Theory 

The  homogeneous  interaction  between  the  beam  and  the  plasma  is 
described  by  the  dispersion  relation 

e(k,oi)  ■  1  +  K  (k,to)  +  Kv(k,w)  =  0  (10a) 

—  P  —  b  — 

where  Kp  and  are  the  longitudinal  dielectric  functions  of  the  plasma 
and  the  beam.  Both  Kp  and  can  be  calculated  for  any  type  of 
distribution  functions  including  collisional  and  finite  size  geometry 
effects  (Briggs  1964).  We  choose  here  models  that  allow  us  to  emphasize 
the  physics  and  avoid  the  mathematical  complexity.  Consistently  with 
Rowland  et  al.  1981  and  Szuszczewlcz  et  al .  1982  we  consider  only  the 
synchronous  Cerenkov  interaction  of  a  slow  beam  wave  with  an  upper 
hybrid  wave  u>Q  of  the  cold  plasma.  In  this  case 
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Wp,  u>k  and  are  the  plasma  frequency,  beam  plasma  frequency, 

2  2  2 

cyclotron  frequency  respectively,  and  k  *  kz  +  k^.  The  poten  *  was 
assumed  to  have  the  form  ^0(^r)  exp[i(kzz-u>t)]  and  thus  the 
propagation  is  axial,  with  k^  determined  by  the  transverse  geometry. 

The  finite  size  beam  reduction  factor  R  enters  through  the  boundary 
conditions  at  the  radii  a  and  b  is  given  in  terms  of  Bessel  functions  by 
(Manickam  et  al.  1975) 


_  a  Yn(klb)  2  2 

R-  2(kia)  b  Y^tJo(kia)  +  Jl(kla)1 


Vkib) 


(11a) 

(lib) 


The  value  tu^R^2  cos9  serves  as  a  reduced  effective  beam  plasma 
frequency.  An  important  aspect  of  the  dispersion  (10)  is  the  abosrption 
of  the  transverse  wavenumber  k^  and  geometry  effects  into  a  single 
parameter  R.  Therefore  from  eq.  (11),  R  is  fixed  when  the  mode  number 
and  the  ratio  —  are  fixed.  Fig.  7  of  Manickam  et  al.  (1975)  shows  the 
values  of  R  as  function  of  for  the  fundamental  mode.  For  large 
Y  »  1,  it  has  a  logarithmic  dependence  approaching  the  value  R  ■  .1. 


In  order  to  determine  the  conitions  for  absolute  instability  in  our 


system  we  follow  the  techniques  developed  by  Bers  (1972),  in  the  weak 


coupling  approximation.  For  the  beam  waves  the  dispersion  relation  is 


2  2 
u),  R  cos  9 

D.(k,w)  =  1 - - - 7=  0  (12) 
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From  this  we  find  the  usual  fast  and  slow  waves  given  by 
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The  slow  wave  is  a  negative  energy  wave  while  the  fast  is  positive,  i.e. 
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where  is  the  wave  energy,  the  beam  wave  amplitude  and  e0  the  free 
space  dielectric  constant.  The  negative  energy  wave  can  couple  in 
synchronous  interaction  with  the  backward  positive  energy  upper  hybrid 
wave  to  produce  an  instability.  Notice  that  for  the  upper  hybrid  wave 

x&k  ) 

D  (k,u)  =■  1  -  2  ■  (13a) 


while  the  wave  energy  is 
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We  examine  now  the  situation  where  the  slow  wave  of  the  beam  interacts 
resonantly  with  the  plasma  wave  u>0.  We  find  a  set  of  coupled  equations 


(14) 
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where  v^,  Vp  are  the  group  velocities  of  the  beam  and  plasma  waves  and 
vb *  vp  are  Phenomenological  damping  coefficients  of  the  two  waves.  uptb 
are  the  usual  normalized  amplitudes  defined  as 
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3D 
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(15) 


(Bers  1972,  Davidson  1970,  Weiland  and  Wilhelrason  1977)  and  the  coupling 


coefficients  C^p,  Cp^  given  by 
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where  Jp^  is  the  perturbed  plasma  current  that  interacts  with  the  beam 
and  Jbp  is  the  perturbed  beam  current  that  interacts  with  the  plasma. 
For  conservative  interactions 
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bp  pb  pb 


(17) 


From  eqs.  (12-16)  we  find  the  growth  rate 
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Absolute  Instability  requires  (Bers  1972) 
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v,v  <  0  (19a) 
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where  L  is  the  system  size  in  the  z-direction.  The  first  condition 
enters  through  the  requirement  that  the  unstable  pulse  encompasses  the 
origin  at  all  times.  The  second  from  the  requirement  that  the  pulse 
growth  exceeds  the  dissipation.  The  last  is  equivalent  to  the  breakdown 
length  Lc  of  an  oscillator  and  implies  that  the  feeback  is  stronger  than 
convective  losses.  Notice  that  in  the  absence  of  wave  reflecting 
boundaries  only  the  upper  hybrid  branch  can  be  absolutely  unstable, 
since  the  lower  hybrid  branch  corresponds  to  a  forward  wave  (i.e.  vbvp  > 

U) 

0 

0).  In  a  plasma  with  —  <  1,  waves  near  the  plasma  frequency  will  be 

e 

convectively  unstable. 

As  mentioned  in  the  introduction  we  associate  the  threshold  of  BPD 

in  the  Johnson  chamber  with  an  absolute  Instability  near  u>e.  For  our 

parameters  the  collislonal  ionization  generates  a  plasma  with 
0) 

>  2,  and  the  collisionality  is  such  that  condition  (19b)  is 
e 

trivially  satisfied.  We  therefore  concentrate  on  eq.  (19c),  for  a 

u) 

0 

backward  wave  in  the  upper  hybrid  range  with  —  »  1.  The  group 

e 

velocity  of  the  waves  vb,  vp  are 
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The  value  of  sln6  can  be  computed  using  (lOe)  and  the  first  root  of  eq. 
(lib),  l.e.  k^b  -  2.4,  giving 
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From  eqs.  (18),  (19c),  (20)  and  (21)  we  find  the  criterion  for  absolute 
Instability  as 
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Iv.  BPP  Ignition  Scaling 

Eq.  (22)  is  the  threshold  condition  for  an  absolute  instability 
near  the  plasma  frequency  in  terms  of  the  plasma  parameters  of  the 
system.  Using  eq.  (9)  we  find  the  current  threshold  condition  as 
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In  order  to  make  further  process  we  have  to  specify  the  value  and 
scaling  of  the  diffusion  coefficient  D.  In  the  low  pressure  regime  it 


was  experimentally  determined  (Szuszczewicz  et  al.  1979)  that  the 
diffusion  coefficient  obeyed  the  Bohm  diffusion  law.  Therefore 
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Using  this  expression  in  eq.  (23)  with  Tg  =>  2.5  eV  we  find 
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For  values  of  g  £  1 ,  the  factor  in  parenthesis  is  between  2-3,  so  that 
to  within  a  factor  of  two  (25a)  reads 
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For  the  standard  parameters  (E,  *  1  keV,  B  =>  1G,  L  *  20  m) ,  we  find  a  at 

D 

the  transition  point  between  the  high  and  low  density  regime 
(pQ  "  15  uT)  the  minimum  current  required  for  BPD  triggering  as 
Ic  -  11  mA.  The  case  analyzed  in  detail  by  Kellog  et  al .  (1982),  i.e. 
p  »  5  uT,  gives  Ic  -  22  mA.  Both  values  compare  favorably  with  the 
observed  values.  While  Bohm  diffusion  is  independent  of  pressure, 
classical  diffusion  in  pressure  dependent,  i.e. 
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We  therefore  expect  that  at  some  value  Dc^  >_ Dg.  From  eqs.  (23)  and 
(26)  we  find  for  Tg  *  2.5  eV 
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For  the  Inflection  point  (i.e.  p  “15  pT)  this  will  give  I  »  14  mA. 

o  c 

Notice,  however,  that  Ic  ~  p^^,  while  rc  ~  "g  rather  than  found 

B 

for  the  low  pressure  regime. 

An  interesting  scaling  results  from  the  above  for  the  dependence  of 
the  minimum  pressure  pQ  for  BPD  as  a  function  of  B.  Namely  pQ  ~  B. 

Since  p^  ■»  15  pT  for  B  =  1G,  we  find  pQ  in  pT  as 

PQ  -  15  B  (28) 

Namely  the  pressure  range  of  the  low  pressure  scaling  Increases  with 

B.  This  is  consistent  not  only  with  Fig.  1,  but  also  with  the  small 

chamber  results  (Konradi  et  al.  1983,  Bernstein  et  al.  1983)  in  which 

B  “  38G  and  the  low  pressure  scaling  was  consistent  with  eq.  27a.  Let 

me  finally  note  that  eqs.  (27a,b)  are  consistent  with  more  recent 

experimental  results  (Lyakhov  et  al.  1982,  Kawashima  et  al.  1982). 

An  approximate  criterion  for  the  ignition  of  BPD  was  given  before 

(Rowland  et  al.  1981,  Papadopoulos  1981)  and  was  found  consistent  with 

the  observations  at  the  J.S.F.C.  tank  (Szuszczewicz  et  al .  1982)  as 
u> 

6 

q—  >  5.  It  is  appropriate  to  comment  on  its  relationship  to  the 
e 

present  more  detailed  considerations.  Referring  to  eq.  (23)  we  note 
that  Ic  ~  sin9  so  that  for  sin9  «  1  Ic  becomes  very  small. 


1  5 


independently  of  other  considerations.  From  eq.  (21)  sin0  «  1 
corresponds  to 


<2,) 

D 


which  Is  the  opposite  limit  from  the  one  considered  before.  Taking 

u>  n  , 

e  2 

b  «>  0(a)  and  using  eq.  (la)  we  recover  the  condition  >  g-  *  5 

*  e  d 

as  an  approximate  condition  at  which  absolute  instability  develops. 

This  criterion  is  extremely  relevant  for  cases  with  preionized  plasma 
such  as  the  ionosphere  at  daylight  conditions  or  high  altitude  (i.e.  F 

Cl) 

0 

peak).  The  —  >  5  criterion  is  a  sufficient  but  not  a  necessary 
e 

condition,  and  accounts  in  a  natural  fashion  for  the  observed  hysterysis 
during  PBD  extinction  (Bernstein  et  al.  1979). 


V.  Summary  and  Conclusions 

We  presented  a  detailed  physical  analysis  of  BPD  threshold  scaling 
based  on  the  conjencture  proposed  by  Papadopoulos  (1981),  that  the  BPD 
threshold  is  associated  with  the  triggering  of  an  absolute  Instability 
near  ojg.  This  conjecture  predicts  a  scaling 

<n 

h  '  pl/2  B‘/2  L 


for  low  pressures  (p  <  p0)  and 


I 
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_3/2  1/2 
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for  high  pressures  (p  >  pQ).  In  deriving  these  scalings  we  assumed  Bohra 
diffusion  for  p  <  pQ  and  classical  diffusion  for  p  >  pQ.  The  minimum 
threshold  current  at  p  =  pQ  is  associated  with  the  transition  from  Bohm 
to  classical  diffusion.  The  predicted  scaling,  the  numerical  values  of 
Ic,  as  well  as  the  dominance  of  Bohm  diffusion  for  p  <  pQ  are  in  good 
agreement  with  the  data  from  JSC  tank  experiment  (Bernstein  et  al.  1979, 
Szuszczewicz  et  al.  1979).  It  is  important  to  restate  some  of  the  key 
assumptions  of  the  theory  and  their  consequences,  since  many  were 
inspired  from  the  JSC  tank  parameters  and  might  not  be  applicable  to 
other  situations. 

(a)  Weak  magnetic  field  in  the  sense  that  at  prebreakdown 

0)  ti) 

—  >_  2.  For  situations  where  —  <  1,  the  toe  waves  lie  in 
e  e 

the  lower  hybrid  branch  of  the  dispersion  curve  which 

corresponds  to  forward  waves  and  therefore  produces 

convective  instability.  This  seems  to  be  the  situation  in 

Boswell  and  Kellog  (1983)  and  the  expected  scaling  should  be 

derived  from  different  considerations.  The  case 

<  we  <  2ng,  requires  special  consideration  due  to  the 

presence  of  strong  cyclotron  damping  which  was  neglected 

here.  Note  also  that  end  plate  reflections  can  produce  an 

absolute  instability  even  in  the  lower  hybrid  branch. 

(b)  Thin  beam  in  the  sense  a  «  b.  For  situations  where 

a  «  b,  axial  losses  could  be  dominant  and  the  more  general 
eqs.  (8)  and  (9)  should  be  used  instead  of  (7). 

(c)  Weak  coupling  limit  eqs.  (14)  are  valid  in  the  limit  where 

the  parameter  «  1. 

e  o 


Otherwise  the  saddle 


point  method  (Le  Queau  et  al.  1980)  should  be  used  to 
determine  the  threshold  condition  given  by  eq.  (22). 

(d)  In  associating  the  threshold  for  absolute  instability  with 
the  BPD  ignition  we  implicitly  assumed  that  the  energy 
deposition  from  the  beam  to  the  plasma  will  produce  ambient 
electron  fluxes  whose  ionization  rate  exceeds  the  beam 
ionization  rate.  This  seems  to  be  clearly  satisfied  in  the 
low  pressure  regime.  However  as  noted  in  Papadopoulos  et  al. 
(1983)  by  increasing  the  presure  there  is  a  limit  at  which 
the  energy  deposition  by  the  beam  plasma  instability  is  not 
sufficient  to  overcome  line  emission  so  that  the  electron 
energy  stays  below  the  ionization  energy.  This  case  will 
have  many  of  the  BPD  signatures  (i.e.  broadening  of  the 
radiated  region)  but  the  electron  plasma  density  will  be 
controlled  by  beam  ionization  only  (i.e.  BPD  without  D) .  The 
small  chamber  experiment  (Bernstein  et  al.  1983)  is  possibly 
indicative  of  such  behavior. 

Before  closing  we  should  comment  on  the  applicability  of  the  above 
concepts  to  electron  beam  injection  in  space.  There  are  two  fundamental 
differences  between  the  laboratory  and  the  ionospheric  beam  injection: 

(i)  The  existense  of  ambient  plasma  with  long  density  gradient 
scales  (i.e.  L  +  °»)  in  the  ionosphere. 

(ii)  The  laboratory  experiments  are  steady  state,  while  the 

vehicle  motion  across  the  magnetic  field  line  can  be  thought 

as  producing  beam  pulses  with  time  length  t  ■  where  U  is 

the  cross-field  motion  of  the  vehicle  (i.e.  1-2  for 

sec 

km 

rockets,  4-8  — —  for  the  shuttle). 
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In  assessing  BPD  for  the  ionospheric  case  we  have  to  ask  the  following 
questions : 

(i)  Is  the  plasma  density  produced  by  the  collisional  ionization 
due  to  the  beam  during  t  «  t,  larger  than  the  ambient 
ionospheric  density. 

(ii)  Is  there  a  beam  instability  at  based  on  the  ambient 
ionospheric  plasma  density. 

(ill)  Is  the  ionization  time  due  to  the  hot  electrons  produced  by 
the  instability  shorter  than  the  injection  time  t  (i.e. 

vion  T  »  1). 

If  the  answer  to  question  (i)  is  yes  and  to  question  (ii)  no,  the 
threshold  condition  is  similar  to  eq.  (23),  with  the  length  L  given  by 
the  plasma  density  gradient  due  to  collisional  ionization  by  the  beam. 

In  view  of  the  short  injection  time,  however,  eq.  (23)  is  sufficient 
only  for  the  beam  plasma  instability  (BPI)  but  not  BPD.  In  order  to 
have  BPD  we  need  in  addition  to  eq.  (23)  a  positive  answer  to  question 

(iii).  The  condition  Vjon  t  »  1  is  equivalent  to  the  Townsend 
condition  (Galeev  et  al.  1976,  Papadopoulos  1981).  If  the  answer  to 
question  (i)  is  no  and  to  question  (ii)  is  yes,  the  BPI  criterion  will 
be  given  by  eq.  (22)  with  nQ  the  value  of  the  ambient  plasma  density. 

BPD  requires  in  addition  vion  T  >  The  above  comments  should  be  only 
taken  as  guidelines.  More  precise  considerations  require  non-linear  BPI 
computations  for  the  evaluation  of  v^on  as  a  function  of  the  beam 
parameters  and  will  be  reported  in  the  future. 
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Figure  1 


Figure  Caption 


Current  threshold  for  BPT)  Ignition  as  a  function  of 
pressure  for  three  magnetic  field  values  (E  =  1.5  kV  and 
L  =  20m). 


